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ABSTRACT 



Most observational techniques in astronomy can be understood as exploiting the various forms of the first-order correlation 
function g^K As however demonstrated by the Narrabri Stellar Intensity Interferometer back in the 1960's by Hanbury Brown 
& Twiss, and which is the first experiment to measure the second-order correlation function g^ 2 \ light can carry more information 
than simply its intensity, spectrum and polarization. Since this experiment, theoretical and laboratory studies of non-classical 
properties of light have become a very active field of research, namely quantum optics. Despite the variety of results in this 
field, astrophysics remained focused essentially on first-order coherence. In this paper, we study the possibility that quantum 
properties of light could be observed in cosmic sources. We provide the basic mathematical ingredients about the first and 
the second order correlation functions, applied to the modern context of astronomical observations. We aim at replacing the 
Hanbury Brown & Twiss experiment in this context, and present two fundamental limitations of an intensity interferometer: 
the requirement of a chaotic light source, and the rapid decreases of the amount of correlated fluctuations with surface's 
temperature. The first of these limitations emphasize paradoxically the fact that the exploitation of g^ 2 ' is richer than what a 
modern intensity interferometer could bring and is particularly interesting for sources of non-thermal light. We also discuss new 
photon-counting avalanche photodiodes currently being developed in Grenoble, and their the impact on limiting magnitudes of 
an intensity interferometer. We conclude by briefly presenting why microquasars in our galaxy and their extragalactic parents 
can represent an excellent first target in the optical/near-infrared where to observe non-thermal light, and test the use of 
in astrophysical sources. 
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1. Introduction 



There is nowadays a revival of interest in the literature 
about the intensity interferometrjf] (LeBohec & Holder 
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|Dravins||2008| |LeBohec eFal 



2008) |de Wit et al.||2008a|b I (see also, for instance, Borra 
2008 |Jain fc Ralston||2008[ on specific applications of in- 
tensity correlations). Intensity interferometry differs from 
the phase interferometry technique (also called amplitude 
interferometry) as it measures correlations between light 
intensities instead of the interferences of electromagneti- 
cal fields. These correlations can also be used to measure 
stellar radii as shown by the Narrabri Stellar Intensity 
Interferometer (hereafter NSII) installed and operated by 
Prof. Hanbury Brown between 1962 and 1972 in Narrabri 



role in astrophysics. It was the first to measure the second- 
order correlation function g^ 2 \ As such it opened the door 
to the observation of phenomena that cannot be explained 
classically, but requiring a full quantum mechanical treat- 
ment. Since then, the so-called Hanbury-Brown & Twiss 
(HBT) effect is used nowadays in various fields of mod- 
ern physics (see for instance |Baym|1998| Alexander|2003 1, 
but not anymore in astrophysics; the main reason being 
its poor sensitivity. 

Hanbury Brown & Twiss explained mathematically 
the equivalence between a classical and a quantum anal- 



in Australia ( Brown|1974[ ) . This experiment has a peculiar 
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1 An interest strong enough for the International 
Astronomical Union to create a dedicated Working Group. 



ysis of the effect they observed (Brown & Twiss 1957 
1958a|b|c l. They emphasized that the observed phe- 
nomenon exemplified the wave rather than the particle 
aspect of light, and claim originally that it has no depen- 
dence on the actual mechanism by which the light energy 
was originally generated. This is not true and this impor- 
tant point is discussed below. 

Technical aspects aside, Hanbury Brown & Twiss suc- 
ceeded at measuring stellar radii with their Intensity 
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Interferometer because of the following fundamental re- 
lationship: 



(i) 



- „(2) 



1 



where g(D and </ 2 ) are the first- and second-order cor- 
relation functions respectively The modulus of g^ cor- 
responds exactly to the definition of the usual visibility 
measured by a Michelson interferometer. The hidden con- 
dition for this equation to be valid is that the light pro- 
duced by the source must be chaotic. Even if this may be 
true for a vast majority of astrophysical sources, it allows 
to draw the boundaries of validity of the technique. 

The possible exploitation of g^ in astrophysics is 
probably far more richer than simply a modern intensity 
interferometer, as shown by the experiments elsewhere in 
physics. Cosmic light sources where there are no inten- 
sity fluctuations, or where the light emission is not (fully) 
chaotic could therefore potentially represent astronomical 
versions of non-classical light sources created in the labo- 
ratory. In this case, Equ.|l]is not valid and g^ will carry 
additional information that is not contained in g^\ This 
paper aims to be a first step to see if meaningful astro- 
physical measurements on such sources can be achieved 
not only to foster progresses on cosmic objects but also 
on light itself. 

In this paper, we review the fundamental meaning of 
the relationships between <?(D and g( 2 \ We aim at giving 
the basic framework to understand why most current tech- 
niques in astrophysics exploit only g^ 1 ', and what means a 
transition towards the use of g^ 2 \ not only through an in- 
tensity interferometer (hereafter II), but also for new and 
so far unaccessible observables. Along the way, the inten- 
sity interferometer of Narrabri can be put in context. We 
discuss the fundamental limitations of an II, and and eval- 
uate them against the capabilities of new avalanche pho- 
todiodes currently being developed in Grenoble. We finish 
by discussing some speculative ideas to go beyond the stel- 
lar II and exploit the second-order correlation function in 
astrophysics. 

2. Mathematical basics 

2.1. The first-order correlation function 

The definition of the normalized first-order correlation 
function reads: 



S (1) (ri,ti,r a ,i 2 ) 



{\E{r x M)\ 2 ){\E(r 2 M)\ 2 ) 



1/2 



(2) 



where E(r,t) is the total complex amplitude of the elec- 
tromagnetical field of the beam light. A measurement of 
g(D evaluates the first-order coherence which is equal to 
the statistical average (denoted by angle brackets) of the 
correlation between electromagnetical fields. Considering 

2 Parenthesis indicate that the number is not an exponent. 



only the temporal coherence which depends on the time 
difference variable r = t% — t\, we have: 



(1) 9 (L> (r) 



(Dm = 



(E*(t)E(t + r)) 
(E*(t)E(t)) 



(3) 



Note that in the case of a plane-parallel wave where only 
the direction parallel to the beam can be considered, r 
might be redefined as r = t 2 — t\ + (z 2 — z{)/c. 

The case of a Michelson interferometer is schematized 
in Fig. [I] (left panel) . Stellar radii are measured through 
the evaluation of the spatial coherence of the light, while 
<i = t% is ensured by delay lines. The evaluation of the 
spatial coherence of the light is a measurement of how 
far two points sources can lie in a place transverse to the 
propagation direction of the light, and still be correlated 
in phase. 

The path difference d\ in Fig. [T] is cancelled out by 
mechanical means to ensure that the evaluation of the 
electromagnetical wave at two different locations remains 
temporally coherent. Somehow, the role of the chaotic light 
is already present in here: if the electromagnetic wave was 
temporally first-order coherent for any time difference t, 
we would not need to equalize the light paths. Another 
possibility to get closer to this situation is to place a high- 
resolution dispersive element which will render each wave- 
length channel closer to a monochromatic source. 

In the case of chaotic light and broad-band observa- 
tions, the temporal coherence time is very short. Hence 
the fringes vanish rapidly with a time difference r = t 2 — t\ 
larger than the light coherence time r c , whatever the de- 
gree of spatial coherence. We can say that our ability to 
see the variation of the spatial coherence when changing 
the distance d between the two telescopes depends on the 
condition that the temporal coherence is fulfilled. 

By definition, the visibility, also called fringe contrast, 
is the modulus of <?(D : 



V 



9 {1 \rxM,r^t 2 ) 



(4) 



In some sense, a (spatial) visibility can be considered 
as the most basic element of an image. It is directly pro- 
portional to the Fourier Transform (FT) of the angular 
distribution of the light intensity on sky (along the pro- 
jected telescopes baseline), according to the van Cittert- 
Zernike theorem (see e.g. Haniff 2007). Although unusual, 
any imaging device can actually be considered as a inter- 
ferometric instrument performing the FT operation and 
transforming the plane wavefront into intensity peaks cor- 
responding to the source positions. Similarly, one can in- 
voke the Wiener-Khintchine theorem to reconstruct the 
spectral density distribution of a source by measuring the 
temporal coherence of the light wave (see below). The in- 
terferences arc cither directly recorded (like in a Fourier 
Transform spectrograph), or a dispersive element (like a 
grating) performs the FT operation and the spectrum is 
directly recorded onto the detector. 

Roughly, one can say that a visibility, a spectrum or 
an image of a celestial source are the different forms of the 
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evaluation of always the same first-order function, g^K In 
the quantum limit, it means measuring the properties of 
single photons (to the contrary of collective properties of 
a photon stream). So far, observational astronomy has not 
made use of g^ (or any higher-order correlation function), 
with the exception of the NSI^] 

2.2. The nature of the difference between temporal 
and spatial coherence 

There is an important difference between the information 
obtained through the measurement of temporal coherence 
compared to that of the spatial coherence. The former al- 
lows to probe the physical processes of light production 
in single or multiple atom systems, while the second is 
related to the spatial " properties" of the emitting system 
itself. As a matter of fact, in the first order, the tempo- 
ral coherence of light is directly related to the spectral 
density of the source, according to the Wiener-Khintchinc 
theorem: 



F(lu) = -ft 

7T 



<7^(t) exp(wjr) dr 



(5) 



where the symbol 3? denotes the real-part of a complex 
number. 

Although being an information about the object it- 
self too, a spectrum provides mostly information on the 
light emission conditions. In other words, the study of light 
production (in any correlation order) is certainly another 
mean to foster progresses in our understanding of cosmic 
objects, but in a manner which is closer to what high- 
energy astrophysics do. When one are interested to light 
production and light properties of higher order or non- 
classical, it is thus natural, because of the fact that Equ[l] 
is a restricting condition, to look for objects where the 
light emission is not fully disturbed by random processes, 
like in star's atmospheres. 

2.3. The second-order correlation function 

By definition, the general form of the (temporal) normal- 
ized second-order correlation function writes: 



(2) _ (E*(t)E*(t + r)E(t + r)E(t)) 



(E*(t)E(t))< 



(6) 



An important overlooked point in the above equation is 
that generally, the complex electromagnetic fields E do 



3 Strictly speaking, the NSII was preceded by two experi- 
ments of the same nature, performed in Jodrell Bank (England) 
under the lead of R. Hanbury Brown. The first one was oper- 
ating in the radio frequencies and was used to measure the 
angular diameter of the Sun, and to study Cassiopeia A and 
Cygnus A. The second one was a crude optical II dedicated 
to prove the feasibility of the concept, by measuring the an- 



gular diameter of Sirius. See Brown ( 1974 1 for the details and 



subsequent references. Historically, the NSII remains however 
the true experiment demonstrating the correlations of intensity 
fluctuations. 



not commute. Therefore, in the general case, the order 
of the terms cannot be rearranged. In the case of chaotic 
light however the cancellation of cross-terms between ran- 
dom relative phases allows this reorganization, and we can 
write: 

(E*(t)E(t)E*{t + T)E(t + T)) 



9 {2 \r) = 



P 



(I(f)I(* + T)) 
P 



(7) 
(8) 



Using Equ. [6] we can show that (Loudon 2000[ 
Equ. 3.7.13, p. 110), for n independent radiative atoms: 

(E*(t)E*(t + T)E(t + T)E(t)) = 

n 

^(Enmit+T^t+^E^t)) + 



E 



+ (E*(t)E*(t + T )Ei(t + r)E {t))} (9) 

Only the terms where the field of a given atom is multi- 
plied by its complex conjugates are kept. All other terms 
vanish because the average of the waves of the different 
atoms with random relative phases are zero. Therefore: 

(E*(t)E*(t + T)E(t + T)E(t)) = 

n(E*(t)E*(t + T )Ei(t + T)Ei{t)) + 

n(n-l){{E*(t)E i (t)) 2 + \(Et(t)W + r))f} 

In the usual case where n» 1, the second term dominates 
largely. We have therefore: 

(E*(t)E*(t + T)E(t + T)E(t)) = 

n 2 {{EtiW®) 2 + \(EUt)E l (t + r))\ 2 } 

Using the definition of g^ and g( 2 \ we obtain Equ[l] 



ff W(r) = l+ 5 (1) (r) 



for n»l 



(10) 



Note that g^ is a real number to the contrary of g^\ 
To recover the phase information, we can still use the 
amplitude itself thanks to the analyticity of the discrete 



Fourier Transform (sec for instance Holmes & Bclcn'kii 



2004 although an image flip degeneracy remains), or the 
high-order correlations (see for instance |Zhilyaev 20081, 



but with the issue of an even lower sensitivity. 

The correlator of an intensity interferometer is de- 
signed to measure the quantity: 

([l(r,t)-l] [l(r + p,t)-l]) 



P 



9 {2) {p)-l 



(11) 



where p indicates the distance vector between the two 
telescopes, as in Fig. [l] and the mean light intensities / 
and I(r,t), are defined as: 



i = (Ht)) 

J(t) = \ce \E{t)\ 2 



(12) 
(13) 



where it is assumed that the light source is stationary (i.e. 
its statistical properties does not change with time). 
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Fig. 1. Working principle of an intensity interferometer compared to that of a Michelson interferometer. On the left 
(Michelson interferometer), a first-order spatial coherence is made, associated with the statistical average of single point 
sources. Delay lines are used to cancel out the path difference d\ and ensure temporal coherence. This cancellation 
must be achieved to a precision that has to be smaller than the coherence length of the light at the given wavelength: 
Ax < c/t c . On the right (intensity interferometer), a second-order spatial coherence is made, associated with the 
statistical average of the correlations between pairs of point sources. Path differences do not need to be cancelled 
to the same precision, since the dominant fraction of this path difference due to azimuthal distance and atmosphere 
fluctuations is identical for both paths. However, the equality of both paths must be achieved to a level corresponding 
to the difference of coherence time between the two Fourier components: Ax < c/(r c i — t c2 ) = - 2 / r n 3> c/r c i 
when t c \ ~ t c2 . This explains the less stringent need in high-precision optical devices and the very poor sensitivity of 
the II against seeing fluctuations. 



2.4. The demonstration by Hanbury Brown of the 
working principle of the intensity interferometer 

In order to clarify what an intensity interferometer mea- 
sures, and how it is related to the above formalism, let 
us briefly review the explanation provided by Hanbury 
Brown (Brown 1974 Chap. 2). The basic principle is 
schematized in Fig. [1 (right panel). Intensity interferome- 
try is a second-order measurement, and instead of evaluat- 
ing the electromagnetic fields in two different locations, it 
evaluates the correlations between pairs of point sources. 
Each points radiate white light and is completely indepen- 
dent of one another. 

A lightwave can be decomposed as a superposition of 
a large number of sinusoidal components, each component 
having a steady amplitude and phase over the period of 
observation but both the amplitude and phase being ran- 
dom with respect to the other components. Let us con- 
sider only one component of this superposition (Hanbury 
Brown calls it a Fourier component) reaching telescope A 
from P\ and another component with different frequency 



reaching the same telescope A from the second point P%- 
We have: 



C\ = Ei sin(u;it + <f>i) 
C 2 = E 2 sin(uj 2 t + 4> 2 ) 



(14) 
(15) 



where E denotes the electromagnetic field amplitude of 
the given component. 

The intensities are transformed into electrical currents 
in the photodetectors. The output current as measured 
behind telescope A is proportional to the intensity of the 
light. Assuming linear polarization: 

I A = K A [Ei sin(wii + 4>i) + E 2 sin(cj 2 t + <fe)] 2 (16) 

where is a constant of the detector. The same Fourier 
component will also reach telescope B 7 with path differ- 
ences g?i and d 2 . Hence: 

I B = K B [Ei sLa(wi(t + dx/c) + <pi) + (17) 
E 2 sin(uj 2 (t + d 2 /c) + 4> 2 )] 2 
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If we expand the above two equations, we obtain each 
time 4 terms. The first two terms are the sum of the light 
intensities from both components, and the third corre- 
sponds to the intensity of the sum of frequencies (uji+uj 2 ). 
All these terms can be filtered out, and we are left with the 
only term of interest here, corresponding to the difference 
of the frequencies — lu 2 ): 

I A = K A E 1 E 2 cos[( Wl - co 2 )t + (4n - cf> 2 )} (18) 
I B = K B EiE 2 cos[(wi - oj 2 )t + (0i - <p 2 ) + 

(widi - io 2 d 2 )/c\ (19) 

These two components are correlated because they 
have the same beat frequency (u>i — w 2 ) and differ only 
in phase by the quantity (u>idi — u> 2 d 2 )/c. If we put for 
simplicity u>\ « uo 2 — to, the correlation, which is the prod- 
uct of I a and I B writes: 



c(d) = K A K B E\El cos -(d 1 - d 2 ) 



(20) 



The phase difference between these correlated compo- 
nents is not the phase difference of the light waves at the 
two detectors but is the difference between the phase dif- 
ferences (or the difference between the relative phases ob- 
served independently at the two telescopes) of the two 
Fourier components. By simple geometry, and because 9 
is a very small angle (sin(#) « 9), we have: 



-,(d) = KaK b E\E\ cos [2ird9/\] 



(21) 



where d is the separation between the two telescopes. 



Finally Brown & Twiss ( 1957 1 have shown that when com- 



puting the total correlation by integrating Equ. [20] over 
all possible pairs of points on the disk of the star, over 
all possible pairs of Fourier components (which lie within 
the optical bandpass) and over all beat frequencies (which 
lie within the bandpass of the electrical filters) we obtain 
that the correlation is proportional to the square of the 
modulus of the visibility. 

Clearly, by no means in this analysis the hypothesis 
of chaotic light is explicitly made. It plays nonetheless an 
important role. 

2.5. The implicit hypothesis of the chaotic nature of 
light in the Narrabri Stellar Intensity 
Interferometer experiment 



As said in the introduction, Brown & Twiss (19571 never 



made the hypothesis that the light source must be produc- 
ing chaotic light. Note that we use here the term "chaotic" 
instead of "thermal", since chaotic light can be produced 
by independent non-thermal individual sources. The sub- 
tle distinction has its importance in the case of radio 
masers, as discussed below. 

It is interesting to note that the year 1957 is also 



the time of the development of lasers (e.g. Schawlow & 
Townes 1958), the properties of which were probably not 
fully known and understood at that time. In their pa- 
per of 1957, Hanbury Brown and Twiss even claim that 



the effect of correlation between intensity fluctuations had 
nothing to do with the mechanism by which the light was 
actually produced. They end their paper by saying: "[...] 
still less does it imply that the photons must have been 
injected coherently into the radiation field. On the con- 
trary, if one wishes to picture the electromagnetic field 
as a stream a photons, one has to imagine that the light 
quanta redistribute themselves over the wavefront, as the 
radiation field, which may be quite incoherent in origin, 
is focused and collimated into beams capable of mutual 
interference; thus the correlation between photons is de- 
termined solely by the energy distribution and coherence 
of the light reaching the photon detectors." 

Interestingly, Hanbury Brown made no reference to the 
nature of light itself when presenting the above analy- 
sis based on Fourier components, noting that this way 
is "freer from conceptual traps". However, in his book 
in 1974 he presented the analysis in a slightly different 
manner, although not emphasizing the radical difference 
it implies for the interpretation of the phenomenon. At 
the beginning of Chap. 3, he did start by stating that 
"white light of the thermal origin has the properties of a 
Gaussian random process" as shown by several classical 
papers. This property is being used further in his analysis 
but never to express how the results could have been dif- 
ferent in the case of non-thermal light. In the formalism 
presented above, it lies precisely at the point where we 
identify Equ. [6] to Equ. [7] 

It is known however today that photons may follow 
different statistics and that it depends precisely on how 
the collections of photons are emitted. This does not in- 
validate the analysis presented above, but it does prevent 
to interpret the correlation seen in g^ as a measurement 
of the visibility |</^|, and hence as an information about 
the angular distribution of light intensity along a base- 
line separating two telescopes via the classical van Cittert- 
Zcrnike theorem. In other words, Equ. [I] is not satisfied for 
non-chaotic light. For instance, it can be shown that, for 
an ideal monochromatic linearly-polarized laser, we have 
g(!)(r) = g( 2 \r) = 1, independent of r. It also means that 
no additional information at all is encoded in higher-order 
functions in chaotic light. In fact, one can show that for 
such type of light only, it is always possible to write the 
order of the correlation function as sums of products 



of the first order function ( Glauber]|2007 p. 115). 



2.6. Photon statistics and non-classical light 

Interestingly, different types of light can be classified 
thanks to <7^(t), or more precisely, depending on their 
value ofg( 2 )(i-~0). 

A typical experimental setup to measure g( 2 \r) is 
sketched in Fig. [2] A counter is started when a photon is 
detected on A, and run until another photon is detected 
on B, which causes the counter to stop. With time accu- 
mulating, an histogram of events can be gradually built 
as function of the duration during which the counter ran. 
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50:50 Beam splitter 



Photons i=> 



O Stop 
O Start 



Counter 



Fig. 2. Sketch of a typical experimental setup to measure 
g( 2 \ Light or photons come from the left and enter the 
50:50 beamsplitter. If a photon is detected on detector A, a 
counter is started. This counter runs until another photon 



is recorded in Detector B (see Fox 2006). It corresponds 
basically to the Narrabri interferometer setup, if the beam 
splitter is replaced by two telescopes recording an intensity 
each. The principle remains however the same. 




1 1.5 
Time delay r/ r c 



2.5 



Fig. 3. Example of </ 2 ' curves for light of coherence time 
r c . In blue dashed a black-body light with gaussian broad- 
ening, in black dotted with Lorentzian broadening. In red 
dot-dashed is shown an example of </ 2 ' curve for anti- 
bunched light whose signature is g^ 2 '(0) < 1. 



By counting the number of events as a function of the 
time difference (r), we measure g^ 2 \r). Typical results 
are illustrated in Fig. [3] In the case of a intensity interfer- 
ometer, the beam splitter is replaced by the spatial extent 
of the light source. 

Quantum optics experiments have shown that g^ 2 \0) 
can be smaller than unity, which is the signature of non- 
classical light. To compare such result to the classical light 
fluctuations, lets adopt a quantum description in terms 
of photons. Intensity fluctuations can be understood as 
"bunches" of photons. In other words, in chaotic light pho- 
tons tends to arrive in groups because the collection of in- 
dividual photon emission is randomly disturbed (by, for in- 



stance, collisions between atoms in the emitting plasma). 
These bunches are characterized by a distribution of ar- 
rival times that is super poissonian, i.e. with a statistical 
variance following: <r 2 (N) — N + N 2 . 

Similarly, in the section of the light beam of an ideal 
monochromatic polarized laser, the arrival time of a pho- 
ton will not be correlated whatsoever with that of another 
one, <?^( T ) remains equal to unity whatever the value of 
t, and the variance is perfectly poissonian: cr 2 (N) = N. 
From a classical point of view, the intensity of such laser 
has no fluctuations at all, and an infinite coherence time. 

The case of g^ 2 '(0) < 1 is often called anti-bunching. 
This cannot be interpreted in a classical sense (the light 
would have a coherence time larger than that of a laser, 
i.e. larger than infinity) and need a quantum description. 
To be exact, anti-bunching and sub-poissonian statistics 
are not exactly the same (see e.g 



Singh 1983 Zou fc 



Mandel||1990l, even if they tend to be satisfied simultane- 



ously. Formally, anti-bunched light is light whose degree 
of second-order coherence satisfies </ 2 ^(t) > g^(0), while 
light with sub- Poissonian statistics satisfies <7^(0) < 1 
(see also LoudonpOOO Chap. 6.5). 

Photon anti-bunching is only one aspect of the quan- 
tum properties of light among many others: squeezed 
states, intricated photons, slow light and so on. It is how- 
ever the first observable that could be used in astrophysics, 
requiring only a large photon collector, and ultra-fast and 
efficient detectors. It is out of the scope of this paper to 
see if other properties can be studied in practice in an as- 
trophysical context. At this stage, we focus on light statis- 
tics and anti-bunching, which still represents truly an un- 
known territory in optical observations of astrophysical 
sources. 

3. The Intensity Interferometer: fundamental 
limitations and new detectors 

It has been presented in other places some considerations 
about using present and future Cerenkov telescope arrays 
to build a modern Intensity Interferometer (see in partic- 



ular |LeBohec fc Holder|2005l |Ofir fc Ribak|2006bl |de Wit 
|et al.||2008b |. It is already known that one of the main 
issues with II is its much lower sensitivity compared to 
an amplitude interferometer, since a second-order effect is 
measured. On the other hand, it is often believed that by 
improving the various parts of the original II experiment 
and multiplying the number of baselines/telescopes, one 
can build a modern interferometer that could open the 
window towards new scientific results. Although it may 
certainly be true in practice, we show below that there 
are at least two fundamental limitations that do not ex- 
ists with a phase interferometer and that must be kept in 
mind when looking for science applications of an II. 

3.1. It works for chaotic light only 

We have exposed above the reasons why an intensity in- 
terferometer ultimately relies on the chaotic nature of the 
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light produced by the source to make measurements of 
stellar radii. This hypothesis of chaotic light is certainly 
verified in the case of thermal light, where the emission 
processes are dominated by random collisions and Doppler 
broadening. Hence, a modern intensity interferometer will 
remain (only) an imaging device of thermal sources. 

Interestingly, second-order coherence has been looked 
for in radio masers, but no deviation from a super- 
poissonian statistics has been found (|Evans et al. 1972 



Moran 1981 1. The only interpretation proposed is that 



even though a maser produce locally coherent light, an 
astrophysical object called "maser" is made of a collection 
of individual coherent sources that have random motions 
relative to each other. 

This example illustrate the obvious fact that the inter- 
pretation of photon statistics in an astrophysical observa- 
tion is not necessarily straightforward and requires some 
a priori knowledge about the source. 

3.2. The hotter, the better 

The second fundamental limitation of an II, to be placed 
on top of the previous one, is that the amount of intensity 
fluctuations, i.e. the "signal" used to measure stellar radii, 
decreases rapidly with the star's surface temperature. In 
other words, the hotter the better. For a thermal source, 
the photon statistics follows the Bose-Einstein distribu- 
tion, whose variance writes: 



(22) 



The first term, called the shot noise, is originating from 
the discrete nature of light (photons) while the second, 
called wave noise, truly represent the fluctuations of the 
energy of the electromagnetic radiation. Hence this latter 
term is the " signal" . However, its contribution for the total 
noise is rather small at optical frequencies. 

Following |Goodman| (|l985| Equs. 9.3-8 and 9.5-19), 



assuming that individual count-products at the output of 
the device for a single counting interval are independent 
from one interval to another, the Signal-to-Noise writes 
for a blackbody: 



V 



III 



(23) 



where V is the first-order visibility, T cxp the exposure time, 
r the shortest integration interval of the detector, and 
8 the degeneracy parameter. The degeneracy parameter 
represents the average number of counts per mode, or per 
single coherence interval of the incident radiation. In other 
words, when (5< 1, shot noise dominates, while 5 ^> 1, the 
contribution of the intensity fluctuations to the statistics 
are dominant. For a blackbody: 

exp(hv/kBT) — 1 ^ 

where T is the source's temperature (see Goodman"| [T985 
Chap. 9.5) . We note that wave noise is by far the dominant 
source of noise in the radio-frequency domain. 
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Fig. 4. Exposure time required to observe a given S/N ra- 
tio with a given visibility and detector integration interval, 
as a function of the stellar temperature, normalized to so- 
lar values. The blue solid curve shows the case for visible 
light (A = 0.55/im), while the red dashed curve is for the 
infrared K-band (A = 2.5/zm). The approximate position 
of spectral types are also shown. 



For a given Signal-to-Noise ratio, visibility and inte- 
gration interval, we have computed the exposure times as 
a function of the star's temperature, and normalized it by 
solar values. The result is shown in Fig. [4] One can clearly 
see that, compared to the observation of a solar analog, 
cooler stars rapidly require much longer exposure times, 
especially at short wavelengths while it is not so critical 
when observing at infrared wavelengths. This limitation 
illustrate the fact that Hanbury Brown & Twiss observed 
mostly hot stars, and it must be taken into account when 
looking for scientific application of a modern II. 



3.3. New photon-counting HgCdTe Avalanche 
Photodiodes and the quantum limit 

The intensity interferometry is the measurement of a 
second-order effect and as such, is poorly sensitive. This 
limitation explains the need for large light collectors, and 
therefore the concomitant interest for this technique with 
the future advent of very large optical telescopes. On the 
other hand, another component not less critical are the de- 
tectors. More efficient and rapid detectors relieve slightly 
the requirements on mirror's size although they increase 
the requirement of isochronicity of the mirrors. As matter 
of fact, they are of central importance when looking for 
light statistics, since the quantum efficiency measures the 
fidelity of this statistics. Using a semi-classical approach 
for photodetection, one can show (e.g. Fox|2006 1 that the 
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observed variance of the statistics of photon arrival times 
writes: 



cr 2 (N) = r] 2 a 2 (n) + 77(1 - 77)7 



(25) 



where the quantum efficiency 77 can be defined as the ratio 
of observed mean photon number N over that of the true 
photon number n reaching the detector. The above equa- 
tion shows that a high quantum efficiency is therefore of 
central importance to have o- 2 (N) as close as possible to 
a\n). 

The other obvious parameter to consider is the small- 
est time interval accessible by the detector, or inversely its 
bandwidth. It can be shown that, for instance in the case 
of chaotic light, the observed photon statistics variance 
writes: 



a 2 (n) 



(7) 



(26) 



where t c is the coherence time of the light, and r the 
exposure time of an individual measurement. If t > r c , 
the variance reduces to a poissonian one. Interestingly, 
the above equation is also valid in the case of the NSII, 
where r c is now the coherence time of the light fluctua- 
tions (which is much longer than the coherence time of the 
thermal light of the observed stars, since r c oc l/(wi— 0J2), 
see above). 

Detectors considered by the time of QuantEYE 
prospective work (a concept of photon-counting instru- 
ment for the E-ELT, see Dravins et al. 2005| and below) 
were limited to the nanosecond accuracy. However, new 
perspectives could be opened by a novel type of polar- 
ized HgCdTe Avalanche Photodiodes developed at the 
Commissariat a l'Energie Atomiquc (CEA) in Grenoble 
( Rothman et al.|2008 ). It appears that these diodes could 
provide for a high gain, a low dark current and a very high 
quantum efficiency (above 95%, J. Rothman, private com- 
munication) a way to detect single photons in the optical- 
near infrared domain with an accuracy of time tagging 
better than 100 picoseconds. 

Fundamentally speaking, one can define an ob- 
servational "photon-counting quantum limit" by using 
Heisenberg's uncertainty principle: 



At AE > h 



(27) 



Using the expected accuracy of the detectors above (80 
picoseconds), and a central wavelength of A ~ 6000 A, the 
limit is achieved with a spectral resolution of R = A/ A A ~ 
40 000. It means that to reach the limit, these detectors 
could in principle be placed behind an echclle grating in 
the optical domain similar to what is installed in UVES, 
the UV- Visual Echelle Spectrograph of ESO's Very Large 
Telescope. 

One can evaluate the impact of such detectors by using 
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Fig. 5. Impact of the new APDs on the limiting magni- 
tude of an II. The comparison is made with LeBohec & 



Holder (2005 1, using a minimal S/N of 3, and a 10 hours 
night integration, an overall reflectivity of R — 0.6, which 
gives a limiting magnitude of about 8 for a telescope area 
of A = 100 m 2 (i.e. d = 11.3 m). The gain with new 
APDs, with an assumed quantum efficiency of 95% and 
a time tagging accuracy of 80 picoseconds is about 2.3 
magnitudes. 

where A is the telescope area in square meters, 77 the pho- 
todetector quantum efficiency, R the overall optics reflec- 
tivity, n the source spectral density in photons per unit 
optical bandwidth, per unit area and per unit time. V 
is the visibility, r the detector and electronic bandwidth, 
and T the total integration time in seconds. For the new 
detectors, we assume a quantum efficiency of 0.95 and an 
inverse bandwidth of 80 picoseconds. The comparison with 
typical detectors expected in Cenrenkov intensity interfer- 
ometers by LeBohec & Holder ( 2005]) with 77 = 0.4 and a 
bandwidth of GHz are shown in Fig. [5] One can see that 
these new APDs would increase the limiting magnitude 
by about 2.3. 

We emphasize that these estimations are somehow op- 
timistic. Brown (1974) observed true S/N systematically 
lower by 25% compared to his theoretical calculations 
given by the above formula. He gave no explanations, but 
identified the most probable cause of this discrepancy as 
the cumulative result of a number of small errors in various 
parameters: overall loss of the optical system, the excess 
noise and atmospheric effects. 



a Signal-to-Noise ratio formula of Brown (1974 see Equ. 4. New scientific questions 



4.30, p50.) taking into account the telescope size and the 
overall quantum efficiency of the system: 



N. 



= ArjRnV 



RMS 



A revival of the intensity interferometry is interesting as 
it may create, among other things, new means to perform 



(28) 



high-resolution angular imaging (de Wit et al. 2008b I. 



However, to our opinion, this regain of interest could be 
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broader, and should encompass the perspectives already 
outlined in the prospective work of the photon-counting 
instrument QuantEYE (|Dravins et al. |[2005 1 proposed for 
the E-ELT (ID 'Odoricol 20051 at that time the ELT was 



called OWL and was expected to have a 100m primary 
mirror), to which the reader is referred for a broad review 
on observational high-speed astrophysics. In that context, 
one must explore the new possibilities offered by true 
photon-counting devices. We note that some experiments 
already started at the Asagio observatory (AquaEYE; 
Barbieri et al.||2007| [Naletto et al.||2007 l. 

With truly new possibilities, we think it is important 
to explore new questions, that do not simply go slightly 
beyond what is accessible by the current (or even foreseen) 
"classical" observational techniques. In this section we at- 
tempt at contributing to these efforts by presenting why 
microquasars in our galaxy, and possibly their extragalac- 
tic equivalent, the quasars, could represent an excellent 
first target in the optical/near-infrared where to observe 
non-thermal light, and test the use of in astrophysical 
sources. 

We note that Hanbury Brown judged perfectly unlikely 
that the coincidence-counting interferometry was of any 
interest for astrophysics. The reason he gives (end of Chap. 
4) is however essentially technical and related to the fact 
that bandwidth must be limited in order to not saturate 
the detector and permit the counting of every single pho- 
ton. The example he took is however that of a zero-th 
magnitude star. 

4.1. Microquasars under the quantum microscope 

Microquasars are short-period X-ray binaries with one of 
its component being a stellar-mass black-hole. These ob- 
jects are the closest relativistic objects to us. They some- 



times show superluminal jets (e.g. Mirabel & Rodriguez 



1994), and produce copious amount of X-rays. Powerful 
and self-collimated jets are produced at the inner re- 



gions of the accretion disk (e.g. Ferreira et al. 2006) 



Microquasars share the same physics as quasars and prob- 
ably gamma-ray bursts (since they share the same phys- 
ical ingredients: a black- hole, an accretion disk and jets); 
they must therefore be understood in the same frame- 
work ( |Mirabel||2004[ ). Microquasars are nowadays studied 
mostly in X-rays because it is where the physics of the 
closest regions around the black-hole can be probed. 

Interestingly, the jets in microquasars are supersonic, 
which means that the particles emitting synchrotron ra- 
diation arc moving faster than the local perturbations, 
possibly avoiding the problem observed in radio masers. 
Moreover, the jets are being produced very close to the 
central black-hole. Therefore, jets could possibly used to 
probe the region around the black-hole similarly to what 
is made nowadays: the accretion disk and its oscillations 
(called Quasi- Periodic Oscillation - QPOs, see for instance 
van der Klis|2005|) can be u sed to estimate black-hole spin 
( |Remillard fc McClintock||2006| ). 




Frequency v (lag 10 Hz) 



Fig. 6. Theoretical Spectral Energy Distributions of a 
microquasar as a function of the frequency v. In red 
the synchrotron radiation of the jet, in blue the multi- 
temperature black-body radiation of the external disk, 
in orange the Bremsstrahlung, in green the Comptonizcd 
synchrotron emission from the inner disk, and in cyan the 



Comptonization of the external disk (see e.g. Foellmi et al. 
2008a|b[ ) + Foellmi et al. 2009). Left: flux (vF v ) expressed 
in flux units . Right: flux expressed in photon count rate. 
Distance of the object is 10 kpc. Optical and NIR regions 
of the spectrum are shown by gray areas. 



The optical/ncar-infrared region of the spectrum is 
well suited for the study of the interface between the jets 
and the accretion disk around the black-hole. Fig. [6] il- 
lustrate a typical theoretical spectral energy distribution 
(SED) of a microquasar with a black-hole mass of 10M Q 



at a distance of 10 kpc (see e.g. Foellmi et al. 2008a I 



(Foellmi et al. 2009, MNRAS, in prep). The left panel 
shows the observed flux of the various radiating compo- 
nents in erg s _1 cm~ 2 . The most central regions of the 
accretion disk around the black hole produce the strong 
X-ray emission. However, the exact same SED expressed in 
terms of photons (right panel) shows a completely different 
situation. In the far-UV and X-ray domain, the number of 
photons is simply much too small to be used for statistical 
studies in a reasonable amount of time. 

Finally, let us mention an example of another quan- 
tum property of the light produced, in theory, nearby a 
black-hole: the Unruh process (a cinematic version of the 
Hawking radiation, Unruh 1976[ ) is expected to produce 
pairs of photons with maximally entangled polarizations 
( |Schutzhold et al.||2006[ ). 



5. Summary and conclusions 

In this research note we have presented the mathemati- 
cal basics about the first and the second order correlation 
functions, and discussed them in the context of astronom- 
ical observations. 

It appears rather clear that detailed physical and quan- 
titative predictions of quantum phenomena in astrophys- 
ical sources accessible with a photon-counting devices are 
still lacking. In the context of the development of large 
Cerenkov arrays and future Extremely Large Telescopes, 
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we think that these questions deserve more dedicated 
studies. 
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